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PREFACE 

The  David  W.  Taylor  Lectures  were  conceived  to  honor  our  founder  in 
recognition  of  his  many  contributions  to  naval  architecture  and  naval  hydro- 
dynamics. Admiral  Taylor  was  a pioneer  in  the  use  of  hydrodynamic  theory 
and  mathematics  for  the  solution  of  naval  problems.  He  established  a tra- 
dition of  applied  scientific  research  at  the  "Model  Basin"  which  has  been 
carefully  nurtured  through  the  decades  and  which  we  treasure  and  maintain 
today.  It  is  in  this  spirit  that  we  have  invited  Prof.  Louis  Landweber  to 
be  a David  W.  Taylor  Lecturer. 

Prof.  Landweber  was  born  in  New  York  in  1912.  He  received  his  Bache- 
lor of  Science  degree  in  Physics  from  the  City  College  of  New  York  in  1932. 
That  year  he  came  to  the  U.S.  Experimental  Model  Basin  at  the  Washington 
Navy  Yard  as  a junior  physicist.  When  the  David  Taylor  Model  Basin  opened 
at  Carderock  in  1940,  Dr.  Landweber  headed  a small  research  group  which 
shortly  afterwards  was  expanded  into  the  Hydrodynamics  Division.  During 
his  years  at  Carderock,  he  continued  to  head  this  division,  and  he  also 
completed  work  on  his  Ph.D.  at  the  University  of  Maryland.  He  was  an  advo- 
cate of  advanced  training  for  the  staff  and  taught  courses  for  the 
University  at  the  Center. 

After  twenty-two  years  of  distinguished  research  at  the  "Model  Basin," 
Dr.  Landweber  left  to  become  a professor  at  the  University  of  Iowa.  He 
continued  his  research  in  ship  hydrodynamics  at  the  school's  Institute  for 
Hydraulic  Research.  During  the  intervening  years  since  leaving  the  Center 
in  1934,  he  has  maintained  close  ties  with  his  colleagues  here  at  the 
Center  and  has  returned  frequently  for  meetings  and  panel  sessions. 


FOREWORD 

After  a lapse  of  twenty-four  years,  it  was  surprisingly  easy  and 
pleasant  to  readjust  for  a month  to  a daily  schedule  of  work  in  a mezzanine 
office  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center. 

I found  Ship  Hydrodynamics  to  be  alive  and  well,  its  problems  being  vigor- 
ously attacked  by  a dedicated  and  talented  staff,  as  in  the  "old  days." 

The  adjustment  was  probably  equally  easy  and  pleasant  for  another 
"alumnus"  and  David  Taylor  Lecturer,  John  Wehausen.  It  is  remarkable  that 
so  many  of  the  alumni  have  stayed  in  the  field  of  Ship  Hydrodynamics,  an 
indication  of  the  challenging  nature  and  attractiveness  of  the  subject. 

To  all  who  visited  me,  to  consult,  gripe,  reminisce  or  simply  to 
educate  me,  who  extended  generous  hospitality  in  accordance  with  the 
fortune-cookie  admonition,  "Take  time  to  play  in  order  to  have  a long 
life,"  who  invited  me  and  arranged  that  day-to-day  living  would  be  con- 
venient and  comfortable  for  my  wife  and  myself,  in  other  words,  to  all 
who  helped  disprove  the  old  adage  that  "one  cannot  go  home,"  our  deepest 
appreciation. 


ON  I RROTAT ION AL  FLOWS  EQUIVALENT  TO  THE 
BOUNDARY  LAYER  AND  WAKE 


INTRODUCTION 

"^Or.e  of  the  simplest  and  most  useful  results  of  boundary-layer  theory 
is  that  the  flow  exterior  to  the  boundary  layer,  which  will  be  assumed  to 
be  irrotational,  is  pushed  outwards  by  an  amount  called  the  displacement 
thickness.  This  suggested,  early  in  the  development  of  boundary-layer 
theory,  that  the  accuracy  of  a boundary-layer  calculation  for  a body  could 
be  improved  by  adding  the  boundary-layer  thickness  to  the  body  dimensions 
and  using  the  predicted  pressure  distribution  on  the  so-thickened  body  in 
a recalculation  of  the  boundary  layer.  The  irrotational  field  about  the 
thickened  body,  including  the  displacement  thickness  of  the  wake,  is  itself 
of  great  interest,  and  numerous  attempts  have  been  made  in  this  manner  to 
calculate  the  effect  of  the  boundary  layer  and  wake  on  the  outer  irrota- 
tional flow. 

The  concept  of  the  thickened  body  gives  an  approximate  model  which  is 
usually  justified  by  its  consistency  with  the  approximations  of  thin 
boundary-layer  theory.  We-'shalt^eview  the  basis  of  this  model  and  suggest 
ways  of  refining  it.  Such  a development  would  be  useful  for  several  current 
problems  of  ship  hydrodynamics,  among  them  the  determination  of  Betz 
sources  in  a method  of  calculating  viscous  drag  by  an  analysis  of  wake 
survey  data,  and  the  investigation  of  the  effect  of  the  boundary  layer  and 
wake  on  wavemaking  resistance. 

We  shall  assume  that  the  boundary  layer  and  wake  (BLW)  are  known, 
and  seek  an  irrotational  model  which  yields  the  actual  outer  flow.  An 
approximate  solution  of  this  problem  has  been  given  by  Preston  [1]  and 
Lighthill  [2]. 

Solution  of  Preston  and  Lighthill  - Two-Dimensional  Case 

Let  us  consider  a two-dimensional  flow  about  a body  of  small  curvature 
in  a uniform  stream  of  velocity  Uq,  on  which  a thin  boundary  layer  is 
present.  Coordinates  parallel  and  normal  to  the  surface  will  be  denoted 
by  (s,n),  and  the  corresponding  velocity  components  by  (u,v).  We  shall 


1 


suppose  that  the  radii  of  curvature  are  so  large  compared  with  the  boundary- 
layer  thickness  6,  that  the  curvilinearity  of  the  coordinate  system  may  be 
ignored.  The  normal  velocity  component  at  n = 6 is  then  given  by 


0 


where  U = u (s,6)  and  6^  is  the  displacement  thickness.  In  the  irrotational 
flow  about  the  body,  with  velocity  components  Uj(s,n),  V^Cs.n),  we  have 
(3Uj/3n)Q  = 0 where  the  subscript  0 denotes  n = 0.  Hence  we  may  put 
UjCs.n)  = Uj.(s),  since 

/3V> 

Uj.(s,n)  =.  U^s.O)  + n y^~J0  =•  U].(s,0)  = ll^s),  0 £ n £ 6 

The  normal  velocity  in  this  irrotational  flow,  at  n = 6,  would  then  be 
given  by  V^(s,6)  =.  Vj(s,0)  + 6(3Vj/3n)g,  and  hence,  since  Vj(s,0)  = 0 and 
3Vj./3n  = - dUj/3s,  we  obtain 


VI(s,6)  =.  - 


(2) 


We  also  assume  that  U = u(s,6)  = U^(s,6)  for  a thin  boundary  layer.  Thus 
the  first  term  of  the  right  member  of  (1)  is  attributable  to  this  irrota- 
tional flow,  and  the  second  term  represents  an  additional  outward  flow  due 
to  the  boundary  layer. 

Put  V'  = (U6,)  and  consider  that  the  "known"  values  of  V'  on  the 

ds  1 

contour  n = 6,  the  edge  of  BLW,  pose  an  exterior  Neumann  boundary-value 
problem  for  determining  a source  distribution  m(s)  on  the  surface  of  the 
body  and  (if  the  body  is  symmetric)  along  the  centerline  of  the  wake. 
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Preston  and  Lighthill  give  the  solution  without  hesitation.  Their  result  is 


m = ^-V(s,6)  = (U6.,) 


2ti  ds 


(3) 


for  which  the  argument  given  by  Lighthill  is:  "This  additional  outflow  is 
exactly  'as  if'  the  irrotational  flow  around  the  body  were  supplemented  by 
the  effect  of  a surface  distribution  of  sources,  whose  strength  ...  is  ..." 

Equation  (3)  does  not  seem  as  obvious  to  the  present  writer.  In  order 
to  derive  it,  one  must  assume  that  the  velocity  U(s)  is  undisturbed  by  the 
distribution  m(s),  an  assumption  which  may  not  be  consistent  with  the 
uniquely  defined  solution  of  the  Neumann  problem.  According  to  this  assump- 
tion, the  body  surface  is  an  equipotential  for  the  distribution  m(s)  (be- 
cause the  tangential  velocity  due  to  m(s)  is  zero) . This  immediately  yields 

_ + 


v'  = 0 


since  V'(s,0_)  = 0 on  the  interior  side  of  the  equipotential  surface  and 
V'(s,0+)  - V'(s,0  ) = 2irm(s) . Then,  from  the  Taylor  expression  V'(s,6)  = 
2^— ) o + ...  =.  V'(s,0+),  we  obtain  the  approximate  solu- 
tion (3)  . 

The  special  case  of  a flat  plate,  (even  in  a nonzero  pressure 
gradient),  however,  confirms  the  approximate  solution  (3)  without  requiring 
the  equipotential  assumption.  Since  the  plate  is  the  limiting  case  of  a 


fv’ 

■J-.  m( 


s) 


m(s) 


Figure  1 

very  thin  body,  the  distributions  m(s)  are  present  on  both  sides.  Applica- 
tion of  the  Gauss  flux  theorem  (the  usual  proof)  im  ‘ediately  gives  Equa- 
tion (3). 
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Another  case  that  yields  (3)  is  that  of  a constant  source  distribution 
on  a circle.  This  gives  a velocity  potential  due  to  a source  of  strength 
2irma  at  the  center  of  the  circle  of  radius  a,  and  hence  V'(s,n)  '=.  V'(s,0)  = 
C2JLmal  = 2TTm.  But  the  local  element  about  the  point  s contributes  7Tm 
to  V' . Hence,  if  we  consider  a nonuniform  distribution  which  is  large  in 
the  neighborhood  of  s,  and  negligible  elsewhere,  we  would  obtain  V'(s,n)  '= 
TTm , half  of  the  Preston-Lighthill  formula.  This  suggests  that,  without  the 
equipotential  assumption.  Equation  (3)  would  give  a good  approximation  only 
for  very  thin  two-dimensional  forms. 

The  distribution  m(x)  on  the  surface  of  the  body  y = f (x) , which 
satisfies  the  boundary  condition  9({>/9ng  = v'(s,0)  can  be  determined  as  the 
solution  of  the  Fredholm  integral  equation  of  the  second  kind,  (see 
Figure  2) 


TTm(s)  + ® m(t) 


1 


in  r _ dt  = V'(s) 
9n  st 


(5) 


in  which  t also  denotes  arc  length  along  the  body  contour.  The  indetermi- 

0 

nacy  of  the  kernel  Jin  rg(_  at  s = t can  be  removed  by  applying,  the 

s 

identity  (at  a smooth  point  of  the  contour) 

j a!;  £n  rst  dt  = j h *st  dt  = 71 
**st 

where  the  complex  vector  r e , extending  from  s to  t,  has  given  the 
pair  of  conjugate  functions 


(6) 
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X-n 


/ 14>  v 

(r  e St)  = X,n  r +i<j) 
'st  > st  yst 


to  which  the  Cauchy-Riemann  Equations  have  been  applied.  We  then  obtain 
from  (5)  and  (6) 


2TTm(s)  + 


n 


3 3 

m(t)  t: — X,n  r -m(s)  tt — X,n  r 

9n  st  3n  st 

s t 


dt  = V' (s) 


(7) 


Since  the  new  integrand  vanishes  at  the  point  s = t,  where  the  largest  con- 
tribution to  the  integral  in  (5)  occurs,  the  form  in  (7)  suggests  that  the 
Preston-Lighthill  formula  should  give  a good  first  approximation. 

From  formula  (3)  for  the  source  distribution,  Lighthill  immediately 
concludes  that  the  displaced  surface,  obtained  by  the  addition  of  the  dis- 
placement thickness,  is  a stream  surface.  This  again  requires  the  assump- 
tion that  all  the  flux  from  the  source  distribution  m(s)  be  outward,  i.e., 
that  the  given  contour  be  an  equipotential  for  the  distribution.  For  then, 
the  Gauss  flux  theorem  gives  the  relation 


where  n(s)  denotes  the  stream  surface  generated  by  m(s).  Comparison  with 
(3)  then  shows  that  n(s)  = 6^. 

A Second  Approximation  for  the  Source  Distribution  - 
Two-Dimensional  Case 

It  appears  to  be  easier  to  derive  a formula  for  a source  distribution 
M(x),  equivalent  to  the  BLW,  for  a thin,  symmetric  body,  when  M(x)  is  dis- 
tributed on  the  axis  of  symmetry,  the  x-axis.  The  given  profile  will  be 
denoted  by  y = f (x) , the  edge  of  the  boundary  layer,  EBLW,  by  y = g(x),  and 
the  profile  thickened  by  the  displacement  thickness  by  y = h(x).  Put 
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df/dx  = tan  y.  For  the  irrotational  flow,  the  velocity  components  in  the 
x-  and  y-directions  will  be  denoted  by  u(x,y)  and  v(x,y). 

We  wish  to  determine  the  source  distribution  M(x)  such  that,  on 
y = g(x),  the  stream  function  ip(x,y)  assumes  the  values 

i p = U(6-61)  ' (8) 


We  have 


(lf)g-  (-v£+u4z^ 


ds  / g 


(9) 


But 


v (x , g)  = v(x, 0)  + g vy(x,0)  = TTM(x)  - g ux<x, 0) 


(10) 


u (x , g)  = u(x, 0)  + g uy(x,0)  = u(x, 0) 


(ID 


since  u^(x,0)  = 0 by  symmetry.  Also,  by  (8),  we  have 


d!)g-s|u(S-{i)!(S) 


(12) 


6 


t 


and  hence,  substituting  (10),  (11),  and  (12)  into  (9)  and  solving  for  M(x), 
we  obtain 

ttM(x)  ' ~ [g (x)  u(x,0)-U(6-61)  ] (13) 

Since,  as  is  seen  from  Figure  3, 

g (x)  = f(x)  + 6 sec  y (14) 

and 

U '=  u(x,g)  sec  Y u(x,0)  sec  y (15) 

we  can  write  (13)  in  the  form 

ttM(x)  = ^ (u(x,0)  [f  (x)+6^  sec  yD  (16) 

Here  we  have  assumed  that  f (x) , 6,  <5^  (u-UQ),  v and  y are  small  quantities 
of  the  first  order  and  terms  of  third  and  higher  order  have  been  neglected. 

For  the  irrotational  flow  about  the  body  without  a boundary  layer,  (16) 
yields  the  well-known,  second-order  approximation  for  a centerplane 
distribution, 


ttM0(x)  = ^ [u(x, 0)  f (x)  ] 


(17) 


Hence,  the  additional  source  strength  due  to  the  boundary  layer  is  given  by 

M(x)  - Mq(x)  '=.  ^ ~ [u(x, 0)  sec  y]  (18) 

The  .result  in  (18)  resembles  the  Preston-Lighthill  formula  (3),  especially 
when  the  body  slope  is  small.  The  distribution  strength  is  approximately 
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doubled  because  a pair  of  distributions  on  the  body  surface  has  coalesced 
on  the  centerplane,  as  was  done  in  the  example  of  the  boundary  layer  on  a 
flat  plate. 

Similarly  to  (14),  the  displacement-thickness  profile,  y = h(x),  may 
be  expressed  by 


h(x)  '=  f (x)  + 6^  sec  y 


(19) 


and  (17)  becomes 

Ml(x)  = i ^ [u(x,0)  h(x)]  (20) 

Comparison  with  the  form  of  (16)  now  shows  that,  to  the  second  order  of 
accuracy,  the  irrotational  flow  about  the  displaced  surface  satisfies  the 
specified  boundary  condition  (8)  at  the  edge  of  the  boundary  layer  and  wake. 

Solutions  for  Axisymmetric  Flow  - First  Approximation 

Let  rg(s)  denote  the  radius  of  a body  of  revolution,  where  s is  arc 
length  along  a meridian  section  of  the  body,  and  r the  radial  distance  to  an 
arbitrary  point.  Let  n denote  distance  from  the  body  along  the  outward 
normal  to  its  surface.  We  shall  also  employ  cylindrical  coordinates,  (x,r) 
where  the  axis  of  symmetry  is  taken  as  the  x-axis,  and  put  rg(s)  = f(x) 
and  tan  y = df/dx.  The  edge  of  the  boundary  layer  is  defined  by  the  surface 


g(x) 


Figure  4 
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r^(s)  = Tq(s)  + 6 cos  y or  g(x)  = f(x')  + 6 (x1)  cos  y (x*)  (21) 


By  a Taylor  expansion  about  x,  with  dy/ds  = - k,  where  k is  the  curvature 
of  the  meridian  profile,  and  with  x’  - x = 6(x’)  cos  y(x'),  g(x)  becomes 


g(x)  = f(x)  + 6 sec  y + 66' 


1 , x2 

tan  y ~ y ko 


tan  y sec  Y + 


(22) 


The  element  of  arc  along  a curve  n = constant  is  h^  ds  = (1+kn)  ds.  The 
equation  of  continuity  is  then 


3(ru) 

3s 


+ [ (1+kn)  rv]  = 0 


(23) 


where  u(s,n)  and  v(s,n)  are  the  actual  velocity  components  in  the  directions 
of  increasing  values  of  s and  n.  We  shall  also  require  the  velocity  com- 
ponents U(s,n)  and  V(s,n)  of  the  "equivalent"  irrotational  flow. 

The  boundary-layer  thickness  of  an  axisymmetric  boundary  layer  is 
usually  defined  as 


6 


1 


r*_ 

f u(s,n) 1 

I r0 

L U(s,6)J 

dn 


(24) 


An  alternative  definition,  in  which  a higher-order  term  is  retained,  is 
given  by 


6* 

I 


■ 


r U(s,6)  dn  = | r [U(s,6)-u(s,n) ] dn  = r^  U 6^ 


Since  r = r^  + n cos  y,  this  yields  the  quadratic  equation 


r0  <S*  + j 6*2  cos  y = rQ  ^ 


(25) 
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and  the  relation 


261  cos  y"| 

1+  — 

L rn  J 


The  equivalent  irrotational  flow  is  determined  by  the  values  of  v(s,6), 
By  (23),  we  have 


(l+k6)  r^  v(s,6)  = 


- / h (ru)  dn  = / h 


[r(U-u)-rU]  dn 


37  - 5 37  ('ou>  - f 37  (u  “s  ^ (27) 


where  U = U(s,6).  In  the  irrotational  flow  about  the  body  without  a bound- 
ary layer,  with  velocity  components  (U^jV^),  (27)  becomes 


0 

(l+k.6)  r^  Vj (s ,6)  = - j [rUI(s,n)]  dn 


“■  - 6 b (roui>  - r k (u! cos  Y) 


since  (9UT/3n)  . = 0,  and  UT(s,n)  = UT(s,0)  + 0(6  ) = UT(s).  If  the 

l n=u  it  i 

boundary  layer  is  thin,  then  Uj(s)  =.  U and  the  first  term  of  the  right 
member  of  Equation  (27)  can  be  interpreted  as  that  contributing  an  addi- 
tional flux  due  to  the  boundary  layer.  We  define  the  additional  normal 
velocity  V'  due  to  the  boundary  layer  by 


(1«S)  t{v  - ^ (VJSj) 
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If  the  body  were  a circular  cylinder  of  constant  radius  rQ,  we  could 
immediately  determine  the  source  strength  m(s)  equivalent  to  the  effect  of 
the  boundary  layer.  Application  of  the  Gauss  flux  theorem  gives 

2nr^V'  = 47T(2nr0m) 


and  hence  by  (28) , 

m = fcH7(u6i>  (29) 

A generalization  of  (29)  can  be  derived  from  the  integral-equation 
formulation  of  the  Neumann  problem  for  the  prescribed  boundary-layer  flux. 
If  P and  Q are  points  on  the  body  surface,  we  have 


2um(P) 


1 


"<Q)^^dvv,<p) 


(30) 


The  singularity  of  the  kernel  when  Q coincides  with  P can  be  removed  by 
means  of  the  relation,  valid  at  smooth  points  of  the  surface. 


(31) 


Thus  (30)  may  be  written  in  the  form 


4TTm(P) 


/ b>^-"<p)^]dv  v,<8-o> 


(32) 


which  suggests  the  first  approximation 
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(33) 


m(P)  = 


1_ 

4tt 


V'(s,0) 


4ur  ds  (rOU6l) 


a generalization  of  (29). 

In  order  to  improve  upon  the  approximate  source  distribution  (33),  we 
shall  now  consider  an  axial  source  distribution  M(x).  For  irrotational 
flow  about  the  given  body,  r = f (x) , the  modified  Munk  formula  of  Landweber 
[3]  gives  the  approximate  solution,  in  terms  of  the  free-stream  velocity  Uq 

M0(x)  = | (!+ki)  uo^f2(x)  (34) 


where  is  the  longitudinal  added-mass  coefficient.  For  the  thickened 
body,  with  r = r^(x)  = f(x)  + 6^  sec  y,  formula  (34)  gives 


Vx)  = 4 (1+kl)  °0  h rl2(x) 


+ \ (l+kx)  U0  ^ (f61  sec  Y + \ 6X2  sec2  y) 


(35) 


If  the  alternative  displacement  thickness  6*  is  used,  the  same  expression 
is  applicable,  with  6^  replaced  by  6*. 

Axlsymmetric  Flow  - Second  Approximation 

We  shall  now  derive  a second-order  approximation  for  an  axial  source 
distribution  for  slender  bodies.  The  given  profile  is  r = tg(x)  = f(x), 
the  edge  of  the  boundary  layer  is  r = g(x),  and  the  displacement-thickness 
profile  is  r = h(x).  The  velocity  components  in  the  x-  and  r-directions 
for  the  equivalent  irrotational  model  are  u(x,r),  v(x,r). 
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Let  us  apply  the  Gauss  flux  theorem  to  the  volume  of  the  body  of 
revolution  generated  by  the  profile  r = f(x),  extending  from  the  nose  to  a 
transverse  section  at  x,  assuming  that  the  flow  about  the  body  is  generated 
by  an  axial  source  distribution  Mq(x).  Since  r = f (x)  is  a stream  surface, 
the  flux  theorem  gives 


4u 


iTTru^x.r)  dr 


(36) 


Here  [uj.(x,r),  Vj(x,r)]  denote  the  velocity  components  of  the  irrotational 
flow  when  the  given  profile  is  a stream  surface,  as  distinct  from  the 
components  (u,v)  of  the  equivalent  irrotational  flow.  Differentiating  (36) 
and  integrating  by  parts,  we  then  obtain 


M0(x) 


1 d_ 
4 dx 


uI(x,f)- 


(37) 
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A 


A 


We  shall  assume  that  f,  6 and  y are  small  of  first  order,  and  neglect 

2 9ux  9(rvj.) 

terms  of  fourth  and  higher  order  in  (37).  Then,  since  r — = r — ~ 

jw  dr  dx 

dM0  3 

= 2r  — — is  small  of  order  f , (37)  yields 


Vx)  ^ [f2uT(x,f)] 


or,  since  Uj(x,f)  = U^s.O)  cos  y =.  Uj(s,6)  cos  y 


MQ(x)  = \ ^ [f2UI(s,6)  cos  y] 


(38) 


Similarly,  the  source  distribution  M^(x)  for  the  irrotational  flow  about 
the  surface  displaced  by  the  displacement  thickness  6*  is  given  by 


Ml(x)  = ^ [ (f+o*  sec  y)2  UI*(s,6)  cos  y] 


where  U^*  denotes  the  s-component  of  the  velocity.  Then,  assuming  that 
Uj(s, 6)  = Uj*(s,6),  we  obtain 


Ml(x)  - MQ(x)  = | ^ [ (2f6*+6*2  sec  y)  UI(s,6)]  (39) 


Similarly,  when  the  surface  of  revolution  is  generated  by  r = g(x)  and 
the  axial  distribution  is  M(x),  the  Gauss  flux  theorem  gives 


4tt  I M(x)  = I 2ng(x)  V(s,6)(l+k6)  dS  + | 2irru(x,r)  dr 
*33  *i)  HO 


or,  by  (27),  and  integration-by-parts  of  the  last  integral. 
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A 


2M<*>  h <U£V  - s S <u£>  lg2u(x,g)j 


The  same  analysis,  applied  to  the  source  distribution  MQ(x),  yields 


2MQ(x)  = 


2 

" 6 h (uif)  ~ f~  57  (ui  cos  Y)  + y 57  [g2Vx’g)] 


and  hence,  assuming  U(s,6)  = U^s.S)  and  u(x,g)  = u].(x,g),  we  obtain 

M(x)  - Mq(x)  = (Uf6x)  (40) 

Comparison  of  (39)  and  (40)'  now  shows  that  M^x)  = M(x)  if 


f6*  + ~ 6*2  sec  y = f6^  (41) 

This  agrees  with  the  definition  of  6*  in  (25)  when  y is  small.  Otherwise, 
the  displaced  surface  should  be  taken  in  accordance  with  (41)  rather  than 
(25). 

A Vorticity  Theorem 

Consider  a steady  mean  flow  of  an  incompressible  fluid  about  a body  at 
rest  in  an  otherwise  unbounded  fluid,  with  mean  vorticity  present  in  the 
boundary  layer  and  wake,  BLW,  bounded  internally  by  the  body  surface  S and 
externally  by  the  surface  T;  see  Figure  6.  The  boundary  conditions  to  be 


Figure  6 
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satisfied  are: 

1.  The  velocity  vector  v on  S+,  the  exterior  side  of  S,  satisfies  the 
nonslip  condition,  v = 0. 

2.  The  flow  exterior  to  T is  irrotational. 

Then  we  have  the  following  theorem: 

The  disturbance  flow  field  exterior  to  S can  be  generated  by  the  vorticity 
in  BLW  alone. 

Proof : 

First  let  us  suppose  that  distributions  of  vorticity  and  irrotational 
singularities  are  present  in  B,  the  interior  of  the  region  bounded  by  S. 

The  velocity  components  normal  to  S induced  by  these  distributions  define 
a Neumann  problem  for  a source  distribution  M on  S,  equivalent  to  the  in- 
terior distributions. 

We  now  have  a source  distribution  M on  S and  vorticity  u)  in  BLW.  The 
tangential  velocity  vt  is  continuous  across  S,  and  hence,  by  the  nonslip 
condition,  v^  = on  the  interior  side  of  S.  Also  the  flow  induced  by  the 
distribi  is  irrotational  within  B.  Hence  v = 0 within  B,  and, 

consequently,  the  _i:mal  component  of  the  velocity  on  the  interior  side  of 
S is  zero.  But,  because  of  the  nonslip  condition,  v^  is  zero  on  the  ex- 
terior side  of  S also,  hence  the  strength  of  the  source  distribution  must 
be  zero.  This  leaves  only  the  vorticity  distribution  to,  as  we  wished  to 
show. 

Verification  with  Stokes  Solution  for  a Sphere 

Stokes  solution  for  a sphere  of  radius  a in  a uniform  stream  in  the 
x-direction,  expressed  in  spherical  coordinates  (R , 0 , 4>) , is  given  by 


Figure  7 
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(44) 


v = U cos  0+2 
K 


vQ  = - U sin  0 + 
0 


(45) 


Here  the  first  and  second  terms  represent  components  due  to  the  uniform 
stream  and  a doublet  at  the  center  of  the  sphere,  and  the  last  terms  are 
due  to  vorticity.  We  shall  now  show  that  the  velocity  at  points  of  the 
,x-axis,  (0  = 0,  R >_  a),  can  be  obtained  from  the  vorticity  outside  the 
sphere,  together  with  the  uniform  stream. 

The  vorticity,  determined  from  the  last  terms  of  (44)  and  (45),  is 
given  by 

2 B 

wR  = 0)g  = 0,  j sin  0 (46) 

R 

At  constant  R and  0,  a vortex  ring  of  unit  strength  induces  a velocity  (by 
the  Biot-Savart  Law),  at  a point  x of  the  x-axis,  given  by 

2 2 „ 

R sin  0 

2 2 3/2 

2(x  -2xR  cos  0+R  ) 


The  velocity  due  to  the  vorticity  in  the  space  exterior  to  the  sphere  is 
then  given  by 


u 


B 


R sin  0 d0  dR 
(x2-2xR  cos  0+R2) 3/2 


(47) 
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With  the  substitution  p = cos  0,  integration  with  respect  to  p yields 


(x2-R2)2  / — — N\+2(x2+R2)(|R-x|-R-x) 

\ | R-x | | R+x | / 

- j ( | R-x | )3-(R+x)3  I dR 


This  is  readily  evaluated  by  considering  separately  the  integrations  from 
a to  x,  and  from  x to  The  former  yields 


u 


1 


A *L 

3 3 


x 


x 

I 


RdR  = - 


Ua 

2x3 


✓ 2 2\ 
(x  -a  ) 


and  the  latter  gives 


OO 


X 


Hence 


ui  + 


3Ua 

2x 


(48) 


in  agreement  with  the  distrubance  velocity  in  (1)  when  0=0,  given  in  (44). 

This  result  indicates  that  the  doublet  at  the  center  of  the  sphere  is 
an  irrotational  equivalent  of  the  negative  of  the  internal  vorticity,  since 
the  sum  of  their  fields  is  zero  exterior  to  the  sphere. 

Some  Applications  of  the  Vorticity  Theorem 

We  have  previously  investigated  irrotational  equivalents  of  the  bound- 
ary layer  and  wake  such  that  the  outer  irrotational  flow  was  preserved.  It 
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was  seen  that  this  could  be  accomplished  in  various  ways,  such  as  a source 
distribution  on  the  given  surface  and  along  the  wake,  a source  distribution 
on  the  edge  of  the  boundary  layer  and  wake,  or  by  the  irrotational  flow 
about  the  displacement-thickness  surface.  In  none  of  these  irrotational 
models  was  it  possible  to  retain  the  given  body  as  a stream  surface. 

Intuitively,  it  appears  desirable  to  match  the  boundary  conditions  on 
both  the  body  and  the  edge  of  the  boundary  layer  and  wake,  in  order  to 
obtain  a more  realistic  irrotational  model.  We  have  seen  that  the  vorticity 
distribution  alone  yields  both  the  nonslip  condition  on  the  body  and  the 
boundary  condition  on  EBLW.  In  the  previous  models  only  the  boundary  condi- 
tion on  EBLW  was  employed.  An  additional  source  distribution,  on  the  body 
or  in  its  interior  is  required  in  order  to  satisfy  the  condition  that  the 
body  remain  a stream  surface.  The  condition  of  zero  tangential  velocity 
would  not  be  satisfied,  but  this  seems  to  be  physically  less  important  in 
an  irrotational  model. 

The  boundary  conditions  on  the  body  surface  S and  on  EBLW  define  a 
Neumann  problem  which  can  be  readily  formulated  as  a pair  or  integral  equa- 
tions to  be  solved  simultaneously  for  a pair  of  source  distributions.  The 
locations  of  the  source  distributions  may  vary,  even  for  a given  body,  as 
has  already  been  illustrated.  If  these  are  taken  to  be  distributions  m(P) 
on  S and  y(Q)  on  T;  one  obtains  the  integral  equations 


2M(P)  - J”(p,)  / 


P PP' 


U(Q')  jj-^dS  , 

9nP  rPQ'  Q 


= - u. 


8x 

0 8n„ 


(49) 


2"U(Q>-  J 


Q QQ' 


m(p,)  7^ — dsp>  = v(Q)  (50) 

d Q P'Q  F 


where  P and  P'  denote  points  on  S,  and  Q and  Q*  on  EBLW,  and  v(Q)  denotes 
the  normal  component  of  the  velocity  at  T,  which  is  assumed  to  be  known. 
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The  Poincare^  Transformation  [4] 

It  is  well  known  that  the  velocity  field  induced  by  a vortex  ring  is 
identical  to  that  of  a doublet  sheet  on  l surface  capping  the  ring.  The 
Poincare/  transformation,  derived  below,  ei ables  more  general  relationships 
between  fields  of  vorticity  and  irrotational  distributions  to  be  obtained. 
Since  the  entire  disturbance  flow  field  is  induced  by  the  vorticity  alone, 
the  velocity  field  of  the  equivalent  irrotational  distributions  would  be 
identical  to  that  induced  by  the  vorticity  in  the  regions  exterior  to  the 
vorticity  domain.  For  the  case  where  the  vorticity  lies  in  the  boundary 
layer  and  wake  of  a body,  the  vorticity  induces  not  only  an  irrotational 
field  exterior  to  its  outer  boundary  T,  but  also  an  irrotational  field 
within  the  body.  Because  of  the  nonslip  condition,  the  induced  velocity 
within  the  body  must  be  zero.  Thus,  in  contrast  with  the  distributions 
previously  considered,  the  Poincare  transformation  offers  the  possibility  of 
matching  the  boundary  conditions  on  both  the  interior  and  exterior  bound- 
aries of  BLW. 

Let  us  suppose  that  vorticity  U)  = curl  v is  present  in  a domain  D, 
bounded  by  a closed  surface  S,  and  denote  the  domain  exterior  to  D by  E. 

Here  v denotes  the  velocity  vector  of  the  fluid  flow.  We  shall  distinguish 
between  a fixed  point  P(x,y,z)  at  which  induced  velocities  are  calculated, 
and  variable  points  of  integration,  Q(£,r]»C).  The  position  vector  from  P 
to  Q is  rpQ  and  has  the  magnitude  rpq- 

The  velocity  induced  by  a vorticity  distribution  can  be  expressed 
either  by  means  of  the  Biot-Savart  Law, 


= i_  f 

4tt  J 3 

rv  /-v 


3 

D ‘PQ 


dT 


or  in  terms  of  the  vector  potnetial. 


v = i-  V x f dT 
VP  4 it  VP  J r 

D ^ 


(51) 


(52) 


20 


We  shall  need  to  distinguish  between  the  vector  operators 


vp  E 1 h + J §7  + R and  v 


1 + j In  + k 


The  Poincare  transformation  is  as  follows: 


X 

f ^ IT  - 

f dS~| 

= VP 

f Pds‘ 

J D P<>  J 

5 J 

r 

JD  PQ 

4 pq  J 

[ 4tt  I _ (P  in  D 
+ < 2tt  > Vp,  < P on  S 
I 0 I P in  E 


where  n is  the  unit  vector  in  the  direction  of  the  outward  normal  to  S. 
First  suppose  that  P is  in  E.  Then  we  have 

) 


f P dS  + f ! 

v*Vp  ~ 

vVp*V  -i- 

JS  D 1 

L PQ 

p p rpqJ 

= V X f dS  + V f v • V — 

P l rPQ  P J P rP0 
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But 


Substituting  (55)  and  (54)  now  yields  (53)  for  P in  E. 

Next  consider  the  case  that  P is  in  D,  and  let  -V-Q  denote  the  volume 
and  SQ  the  surface  of  a small  sphere  of  radius  rQ  about  P.  Then  we  have 


V 


P 


X 


(56) 


where  D'  - D - VQ.  The  last  integral  in  (6)  is  proportional  to  the 
velocity  induced  by  the  vorticity  within  according  to  (52)  and  hence 
must  vanish  as  the  radius  of  the  sphere  approaches  zero.  Also,  by 
Equations  (54)  and  (55) , we  have 
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and 


The  convention  on  the  positive  sense  of  n requires  that  it  be  directed 
inwards  on  the  sphere  S^.  For  the  term  in  (57),  we  have 


b0  S0 


(nxv)  2 „ 

~r~  ro 


so  0 


*•/ 55 


• vdfi  + Attv, 


(59) 


where  dS  = r^  dfi;  and  in  (58), 


VD  f ^dS  = - f n n • 

' -t  rPQ  1 


vdH 


(60) 


Also  we  have 


lim  V 

V° 


/ r!dT 

(61) 

V PQ  JD  PQ 

since  the  velocity  field  of  a volume  distribution  of  sources  is  continuous. 
Collecting  the  results  in  Equations  (56)  to  (61)  now  yields  (53)  for  P in  D. 
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Finally,  when  P is  a "smooth"  point  of  S,  i.e.,  a point  at  which  the 
tangent  to  S is  continuous,  we  introduce  a small  hemisphere  about  P,  of 
radius  r^,  and  apply  the  Poincare  transformation  to  the  so-diminished 
volume  D' , and  to  the  bounding  surface,  consisting  of  the  hemispherical 
surface  and  the  remainder  S'  of  S.  The  proof  is  similar  to  that  for 


P 


P in  D.  In  Equations  (57)  and  (58),  we  need  only  to  replace  S by  S',  and 
in  Equation  (59)  4irvp  becomes  2iTvpr  or,  indeed,  avp  if  P is  a corner  point 
of  S of  solid  angle  a.  Instead  of  (61) , we  need 


lim  V x f pdS'V  x f 7^  dS 
rn-0  Jc(  PQ  Js  PQ 


and 


u»  Vp  f PdS-Vp  fp 

rn+0  I,  PQ  1 PQ 


dS 


which  are  verified  by  observing  that  the  velocity  fields  of  the  surface 
distributions  of  vorticity,  n * v,  and  of  sources,  v • n,  are  continuous. 


Equivalent  Irrotat tonal  Flow  From  Poincare  Transformation 

Let  us  apply  the  Poincare  transformation  (53)  to  a case  in  which  the 
flow  is  entirely  due  to  the  vorticity  in  a domain  D,  bounded  by  a closed 
surface  S.  We  shall  seek  to  express  the  velocity  at  a point  P of  E as  the 
gradient  of  a velocity  potential. 
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The  first  term  of  (53)  is  seen  to  give  the  velocity  vvp  due  to  vortic- 
ity  in  the  form  of  the  vector  potential  (52),  as  4iTVyp.  Since  the  assumed 
flow  is  solenoidal,  the  third  term  of  (53)  is  zero.  We  then  have,  from 
(53), 


4i tv  = V x 

VP  P 


(62) 


This  expresses  the  velocity  in  terms  of  that  induced  by  a source  distribu- 
tion of  strength  v^Mtt  and  by  a vortex  sheet  of  strength  n x v = vgt,  both 
on  S.  Here  we  have  expressed  the  velocity  vector  v = nv^  + svs>  where  s 

is  the  unit  vector  in  the  direction  of  the  projection  of  v on  the  tangent 
plant  at  Q.  Then,  putting  n x s = t,  we  obtain  the  form  given  above. 

In  order  to  express  the  field  of  the  vortex  sheet  as  the  gradient  of 
a potential,  let  us  define  a function  <|>^ , harmonic  in  D,  which  on  S 
satisfies 


n x = n x v,  v^  = V<t>^ 


(63) 


We  have  then 


35"  = Vls  = Vs  °n  S (64) 


and  hence,  by  integration  along  s,  the  values  of  on  S may  be  presumed  to 
be  given  in  terms  of  the  known  values  of  vg.  Thus  the  boundary  condition 
(64)  sets  up  a Dirichlet  problem  for  determining  cfi^. 

Since  V x = 0,  application  of  the  Poincare  transformation  to  v^, 
taking  (63)  into  account,  yields 


(65) 
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Hence  (62)  becomes 


VP 


— V 

4ir  P 


f",(V 

i rPQ 


•v) 


dS 


(66) 


a velocity  field  due  to  a source  distribution  of  strength 

a = n • (v-v^Mir  (67) 

The  form  of  a shows  that  n • is  the  normal  component  of  the  velocity  on 
the  interior  side  of  the  source  distribution  on  S.  Uniqueness  of  solutions 
of  Neumann  problems  on  S then  shows  that  v^  is  the  irrotational  velocity 
field  in  D associated  with  the  source  distribution  0 on  S. 

An  alternative  source  distribution  on  S can  also  be  found  directly  as 
the  solution  of  the  exterior  Neumann  problem  for  the  given  values  of  v • n 
on  S,  as  was  done  previously  in  considering  the  displacement  effects  of  the 
boundary  layer.  Applied  to  the  vorticity  field  BLW  in  the  flow  about  a 
body,  the  present  approach  requires  that  a composite  bounding  surface, 

S = S + T + A 
J5 

be  used,  where  S_  is  the  surface  of  the  given  body,  and  A is  the  transverse 

D 

surface  of  the  truncated  wake.  If  the  transverse  section  is  taken  suffi- 
ciently far  downstream,  the  effects  of  the  source  distribution  on  A may  be 
neglected,  and  the  Neumann  problem  could  be  formulated  as  a pair  of  simul- 
taneous Fredholm  integral  equations  of  the  second  kind.  The  resulting 
source  distribution,  however,  would  not  coincide  with  that  given  in  (67). 

In  the  present  treatment,  the  value  v • n = 0 on  S is  preserved,  whereas, 

D 

in  deriving  (67),  the  tangential  component  v x n = 0 on  S , in  accordance 

D 

with  the  nonslip  condition,  was  preserved.  In  the  former  case,  S remains 

B 

a stream  surface,  in  the  latter,  an  equipotential. 


Betz  Method  for  Determining  Viscous  Drag 

Betz  [5]  and  Landweber  and  Wu  [6]  use  equivalent  irrotational  flows  to 
derive  formulas  for  the  viscous  drag  of  a body  in  terms  of  measured  values 
of  pressure  and  velocity  at  a transverse  section  of  the  wake.  A refine- 
ment of  these  derivations,  in  which  additional  wake  characteristics  are 
taken  into  account,  will  now  be  presented. 

The  body  is  taken  at  the  center  of  a circular  channel  of  large  radius, 
and  is  at  rest  in  a uniform  stream  of  velocity  U in  the  positive  x- 
direction.  The  disturbance  velocity  components  in  a rectangular  (x,y,z) 
coordinate  system  are  (u,v,w),  and  p denotes  the  pressure. 


We  select  a control  surface  consisting  of  the  transverse  sections  AB, 
far  ahead  of  the  body,  CD  or  S a moderate  distance  behind  it,  and  the 
portion  of  the  channel  wall  lying  between  these  sections.  On  the  section 
AT,  designated  S^,  the  pressure  is  the  constant  p^,  and  the  velocity  is 
(Uq,  0,  0).  Application  of  the  momentum  theorem  to  this  control  surface 
yields  the  expression  for  the  body  drag  D, 

D = f {pQ-p-p  [(U0+u)2-UQ2]}  dS  (68) 

in  which  p is  the  mass  density  of  the  fluid.  If  the  wake  is  turbulent, 
Reynolds  stresses  will  be  present,  but  these  can  be  taken  into  account  most 
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efficiently  by  averaging  the  resulting  expression  derived  for  the  drag. 
In  terms  of  the  total  heads,  defined  by 


Pg  Hq  = pQ  + j p UQ2,  pg  H = p + | p [ (U0+u)2+v2+w2]  (69) 


where  g is  the  acceleration  of  gravity,  (68)  becomes 


D'll 


pg  (Hq-H)-  j p [(U0+u)2-Uq2-v2-w2]  } dS 


(70) 


We  now  consider  an  equivalent  irrotational  velocity  field  (Uq  + u^.v^, 
w^),  with  pressure  p^,  generated  by  a volume  distribution  of  sources  of 
strength  u in  BLW,  such  that  (u^.v^.w^)  = (u,v,w)  on  T,  the  outer  boundary 
of  BLW.  We  again  apply  the  momentum  theorem,  to  the  flow  within  the  same 
control  surface  generated  by  this  distribution  of  sd-called  Betz  sources, 
to  obtain  the  expression  for  the  force  on  the  sources  within  the  control 
volume , 


D 

s 


J t(Vui)2-uo2-vi2"wi2l  ds 

•'s 


(71) 


the  term  corresponding  to  the  difference  in  heads  vanishing  since  the 
field  is  irrotational. 

Another  expression  for  the  force  on  the  Betz  sources  is  given  by  the 
Lagally  formula 

Dg  = - 4np  J ^(Uq+Uj^)  dT  (72) 

-V- 


28 


where  Af  denotes  the  control  volume.  This  is  not  the  usual  application  of 
the  Lagally  theorem,  which  gives  the  force  on  a closed  body.  Here  it 
represents  the  reaction  on  the  Betz  sources  due  to  the  flux  of  momentum  of 
the  source-generated  discharge. 

Since  the  flux  through  the  bounding  surface  T is  the  same  with  Betz 
sources  as  in  the  solenoidal  (divergenceless)  real  flow,  the  difference  in 
flux  for  the  flows  through  the  area  A of  S intersected  by  the  wake  is  at- 
tributable to  the  Betz  sources.  The  Gauss  flux  theorem  then  yields  the 
formula 


/ <"r“> 


ds 


(73) 


Far  downstream,  the  section  S will  be  denoted  by  S and  the  wake  area 

J CO 

by  A^.  In  terms  of  the  flux  Q across  this  wake  area,  we  have  the  well- 
known  formula  for  the  drag. 


D = p UqQ,  Q 


(74) 


A similar  expression  for  the  force  on  the  Betz  sources,  obtained  from  the 
asymptotic  form  of  (71),  is 


= - P Ur 


dS 


But 


j Ul  dS  ■ J 


(u^-u)  dS  + 


/■ 


dS 


and  since  u^  = u except  within  the  wake  and,  by  continuity, 
we  obtain 
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(75) 


D = - p 

S°° 


uo  J <V" 


) dS 


and  hence,  from  the  asymptotic  form  of  (73), 


U°I 


Dgoo  = - 4ttP  Uq  | n dT 

BLW 


(76) 


Comparison  with  the  asymptotic  form  of  (72)  now  yields 


/ ^ U1 
BLW 


dT  = 0 


(77) 


From  the  asymptotic  form  of  (73),  we  can  also  show  that 


I 


4tt  I p dx  = Q 
BLW 


(78) 


which  implies,  by  (74)  and  (76),  that 


D = - D 

goo 


(79) 


This  is  derived  by  neglecting  in  comparison  with  u in  (6) , since  the 
Betz  sources  are  concentrated  near  the  body,  so  that  diminishes  as  the 
inverse  square  of  the  distance,  while,  for  a turbulent  wake,  |u|  diminishes 
as  the  inverse  2/3-power  of  the  distance. 

Put  V + V = V and  define  a mean  value  u^  of  u^  by 

p dx  = I p u1  dx  (80) 

v 
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We  then  obtain,  from  (77), 


f P dx  = J p u1  dx  = J 

Jxx  XI  » «.  I 


p dx 


But,  by  (73)  and  (78), 


M dx  = M dx  _ M dx  = a_  _ 1_  f (Ui_u)  ds 

V BLW  W *A 


Then,  by  (74), 


J M U1  dT  = 4?  " J (U1"U)  dS 


(81) 


Hence,  by  (73)  and  (81),  (72)  becomes 


D = 
s 


J (uru)  (U0-Kl1 


) dS  + — D 
0 


(82) 


A formula  for  the  viscous  drag  D can  now  be  obtained  by  subtracting 
the  expression  for  Dg  in  (71)  from  that  for  D in  (70)  and  then  substituting 
for  Dg  from  (82).  Observing  that  the  resulting  integrand  is  nonzero  only 
over  the  wake  area  A,  we  obtain  the  result 


i- 


- / 

ui  i 


[2g(HQ-H)+(u1-u) (u1+u-2u1)+v2+w2-v12-w12]  dS  (83) 
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In  applying  (83),  total  head  tubes  which  measure  a total  head  H , 


given  by 


Pg  Hra  = p + f P [ (U0+u)2+A(v2+w2) ] , 0 < A < 1 


(84) 


where  A is  a calibration  constant,  (see  Reference  [9]),  and  velocities 


1/2 


u = - (H  -p)  - U- 

m L p m J 0 


(85) 


are  used  for  H and  u.  In  terms  of  H and  u , with  the  small  difference 

m m 

between  u^  and  u neglected  in  higher-order  terms,  (83)  becomes 


d = -£^-  f 

1-u, /U_  J 


1/U0  ~A 


[2g(H  -H  )+(u  -u  ) (u ,+u  -2u  ) 
urn  lm  im  1 


2 2 2 2 
-v^  -w^  +A (v  4w  ) ] dS 


(86) 


When  the  wake  is  turbulent,  the  mean  value  of  D,  obtained  by  replacing 
(u,v,w)  by  (u+u' ,v+v ' ,w+w' ) and  averaging,  where  (u',v’,w')  denote  the  com- 
ponents of  the  turbulent  velocity  fluctuations,  becomes 


D=  — ^ — j |"2g(H  -H  ) + (u  -u  ) (u  +U  -2u  )-v  2-w  ^ 
« / I U m X in  X m XXI 

I"ul/U0  A 


+A  (v2+w2)  -u,2+A  (v'2+w'2)  dS  (87) 


where  H and  u now  refer  to  their  mean  values.  The  Reynolds  stress  terms 
mm  J 

2 

combine  into  (2A-l)u'  for  isotropic  turbulence,  and  would  hence  be 
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negligible  for  X = 1/2.  Jin  Wu's  measurements  [7]  indicate  that  the  tur- 
bulent stresses  would  contribute  about  2 percent  to  the  calculated  drag 
with  X = 0;  but  his  actual  value  was  X = 0.5,  with  which  the  turbulence 

terms  became  negligible.  We  shall  assume  that  the  turbulence  stresses  in 

2 2 

(87)  may  be  neglected.  We  shall  also  neglect  the  terms  - v.  - w.  + X 

2 2 11 
(v  +w  ) since  these  are  small  and  partly  cancel  each  other. 

Still  unknown  are  the  terms  and  u^  in  (87).  For  estimating  y^, 

we  shall  assume  that  u.  depends  upon  z alone,  «.nd  is  given  by  u.  = u (z), 

1 1 hj 

the  measured  value  of  u at  the  wake  boundary  T.  According  to  the  defini- 
tion of  u^  as  a mean  of  u^  in  V , we  observe  that  the  mean  is  weighted  by 

the  value  of  the  source  strength  y which,  together  with  u^,  diminishes  to 

zero  as  x -*■  °°.  This  suggests  that  the  values  u,  = 0 and  u.  = u at  A 

1 I t 

can  be  used  to  obtain  bounds  for  the  drag  formula,  the  "true"  value  lying 
closer  to  the  bound  given  by  u..  = u . Because  u.  occurs  both  in  the 
integrand  and  in  the  denominator  of  the  expression  for  the  drag,  it  is  not 
immediately  evident  which  of  the  two  bounds  is  the  larger.  Denoting  these 
bounds  by  when  u^  = 0 and  when  u^  = u^,,  and  applying  the  afore- 
mentioned approximations,  we  obtain  from  (87), 

D1  ' 2 I dS  <88> 

and,  with  u_  + u - 2u„  replaced  by  u - u„  in  the  last  term, 

E m E 1 m E 


The  latter  form  was  given  by  Tzou  and  Landweber  [8]. 

That  < D2  is  indicated  by  the  following  argument.  Since  the  source 
strength  y represents  the  displacement  effect  of  the  boundary  layer  and 
wake,  it  is  a positive  quantity  of  total  strength  given  by  (78).  Conse- 
quently, according  to  (77),  u^  cannot  be  of  one  sign  throughout  BLW.  Over 
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most  of  the  boundary-layer  edge  of  T,  u_  > 0.  In  the  region  of  the  rear  of 
the  body  and  the  near  wake,  there  would  be  some  pressure  recovery  (especial- 
ly if  separation  has  not  occurred)  and  u_  on  T would  there  become  negative, 
and  gradually  approach  zero  with  increasing  downstream  distance.  This 

suggests  that  u in  (89)  and  u.  in  (86)  are  negative.  Since  the  difference 
E J- 

in  heads,  H„  - H , contributes  about  90  percent  of  the  magnitude  of  the 
U m 

integrand  in  (88)  and  (89),  and  u <<  U..,  (89)  may  be  written  as 

E U 


°2  ' f (1+  / 2*<W  dS  - 2 | <V“: 


)2  dS 


which  yields 


°2  - °1  N P J ® (H0~Hm)  dS  " P J UE(UE-Ur, 


) dS 


• r 


i*(ho-h«)-uo(ue-u«)i  dS 


(90) 


Applying  the  expression  for  in  (85;  and 


Pg  H0  = PE  + J P(W2 


2 2 

in  (90)  and  neglecting  the  term  u - u , we  obtain 

E m 


“2  * “1  ■ % l ^ 


dS  > 0 


(91) 


34 


since  u < 0 at  A and,  according  to  Jin  Wu's  data  [10],  p is  a maximum  at 

in 

the  center  of  the  wake,  so  that  p - p > 0.  The  data  given  in  Reference 
[8]  indicate,  however,  that,  at  a section  at  0.6  of  the  length  behind  the 
stern  of  a ship  model,  u^,  is  zero  within  the  accuracy  of  the  measurements. 
This  indicates  that  the  simpler  expression  (88)  is  suitable  for  computing 
the  viscous  drag  from  wake-survey  data. 
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